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Direct determination of the neutrino mass is at the present time one of the most important aims 
of experimental and theoretical research in nuclear and particle physics. A possible way of detection 
is through neutrinoless double electron capture, OuECEC. This process can only occur when the 
energy of the initial state matches precisely that of the final state. We present here a calculation 
of prefactors (PF) and nuclear matrix elements (NME) within the framework of the microscopic 
interacting boson model (IBM-2) for 164^^.^ iso^ pp 

calculate the expected half-lives and obtain results that are of the same order as those of fT^ 
decay, but considerably longer than those of Qvj3~I3~ decay. 

PACS numbers: 23.40.Hc,21.60.Fw,27.50.+e,27.60.+j 


I. INTRODUCTION 


The question whether or not the neutrino is a Majo- 
rana particle and, if so, what is its average mass remains 
one the of the most fundamental problems in physics 
today. In a previous series of papers we have calcu¬ 
lated phase-space factors (PSF) and nuclear matrix el¬ 
ements (NME) for 2i'f3~f3~ processes [H-Q, 

and for 0 1 /EC0'^, and 2vEC0^ and 

2vECEC processes [ij. The neutrinoless electron 
capture OvECEC was not calculated due to the fact that, 
in general, one cannot conserve energy and momentum 
in the process 


{A,Z)+2e- ^ {A,Z -2). (1) 

However, conservation of energy and momentum can oc¬ 
cur in the special case in which the energy of the ini¬ 
tial state matches precisely the energy of the final state. 
This situation was first discussed by Winter 113 and 
subsequently elaborated by Bernabeu et al. [1^, who 
provided estimates of the inverse life-times with simple 
NME and PF. The work of [I^ stimulated many exper¬ 
imental searches and additional theoretical papers |l4l - 
[^ . The process was also termed resonant neutrinoless 
double electron capture, ROvECEC. The matching con¬ 
dition can occur either for g.s. to g.s. transitions or for 
transitions between the g.s. and an excited state in the 
final nucleus. The precise matching condition is an ex¬ 
ceptional circumstance which may or may not occur in 
practice. A slightly less stringent condition is that the 
decay occurs through the tail of the width of the atomic 
initial state as shown schematically in Fig. [T] For this 
process, depicted in Fig. [21 the inverse half-life can be to 



FIG. 1. (Color online) Schematic representation of neutrino¬ 
less double electron capture. 


a good approximation factorized as 




where is a prefactor depending on the probability 

that a bound electron is found at the nucleus (see Sect. 
II), nuclear matrix element, and /(m^, Uei) 

contains physics beyond the standard model through the 
masses rrn and mixing matrix elements Uei of neutrino 
species. For light neutrino exchange 
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FIG. 2. Neutrinoless double electron capture where v = v 
can be either light or heavy neutrino. 
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while for heavy neutrino exchange 
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The last factor, often written as 
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is the figure of merit for this process. Here A = 
\Q — B 2 h — E\ is called the degeneracy parameter, where 
r = Tej + Tf .2 is the two-hole width and B 2 h is the 
energy of the double-electron hole in the atomic shell 
of the daughter nuclide including binding energies and 
Coulomb interaction energy. The importance of the 
Coulomb interaction energy of the two holes was investi¬ 
gated in Ref. M- Obviously the maximum value of F is 

pmax _ 

Since T is of the order of eV one needs to find nu¬ 
clei or nuclear states such that A has the smallest pos¬ 
sible value. This requires an accurate measurement of 
the Q-value. Recent imp rovements in measurements of 
mass differences [3, [H iiS, m [p have ruled out 
many of the initial candidates. In this paper, we report 
calculations of the PF and NME of five remaining can¬ 
didates, 124Xe-^124q.g*^152Qjj_^152gjj^^ 156j3y^l56Qd* ^ 

i64Er_^i64p)y^ and isoW^'isopf, where the star denotes 
excited state. (In addition to these, other candidates re¬ 
main but with spin and parity of the final excited state 
unknown. Theses cases are not discussed here.) The en¬ 
ergetics for the 5 cases of interest are shown in Table ID 
The case ^^^Gd—decay is also illustrated in Fig. 

El 


FIG. 3. Energetics of ^®^Gd—double electron capture. 


II. PREFACTORS 


In the calculation of the prefactor, PF= Gn,9^'^, we 
follow the theory described in our previous papers [lOj , 
in particular that of 0 for electron capture (EC). The 
captured electrons are described by positive energy Dirac 
central field bound state wave functions, 


(l*) 




( 6 ) 


where n' denotes the radial quantum number and the 
quantum number k is related to the total angular mo¬ 
mentum, Jk = |k| — 1/2. The bound state wave functions 
are normalized in the usual way 



(r)dr = 





dr = 1. 


( 7 ) 


They are calculated numerically by solving the Dirac 
equation with finite nucleon size and electron screening 
in the Thomas-Fermi approximation 

For the calculation of electron capture processes the 
crucial quantity is the probability that an electron is 
found at the nucleus. This can be expressed in terms 
of the dimensionless quantity 


^ 

4^(mec2)3 




( 8 ) 


{gn'JR))^ + {fu'AR)y 


where oq is the Bohr radius oq = 0.529 x 10“®cm and we 
use for the nuclear radius R = 1.2A^/^fm. For capture 
from the iX-shell n' = 0, k = —1, lS' 1/2 while for capture 
from the L/-shell n' = 1, k = — 1, 2Si/2- The PF is given 
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TABLE 1. Double electron captures considered, the Q-value of the decay, energy of the resonant state in the daughter nucleus, 
capture shells, energy of the double-electron hole, two-hole width, and resonance enhancement factor. 


Decay 


(5-value(keV) 

F(keV) 

Q - F(keV) 

Shells 

i?2h(keV) 

A(keV) 

r(keV)^ 

(mec2)F 

124-v^ 

54 ^^70 

V 124rp =(= 
^52 -*-^72 

2856.73 ±0.12^ 

2790.41 ±0.09 

66.32 

K-K 

64.457^ 

1.86 

0.0198 

2.92 

IfGdss 


55.70 ±0.18^ 

0.0 

55.70 

K -Li 

54.795^ 

0.91 

0.023 

14.38 

efDygo 

^64®Gdg2 

2005.95 ±0.1^ 

1988.5 ±0.2 

17.45 

Li — Li 

16.914^ 

0.54 

0.0076 

13.52 


■^66^^798 

25.07 ±0.12® 

0.0 

25.07 

Li-Li 

18.259® 

6.81 

0.0086 

0.095 

^fWioe 

^^i°Hfi08 

143.20 ±0.2? 

0.0 

143.20 

K-K 

131.96£1 

11.24 

0.072 

0.29 


Ref. 
Ref. 
Ref. 
Ref. 
<= Ref. 
f Ref. 
s Ref. 
h Ref. 




TABLE II. Prefactors for 

neutrinoless double electron cap- 

ture. 



Decay 




This work 

Ref. 

ifXero 

2.57 


gl^Gdgg ^-g^^Smgo 

1.46 

1.67 

li^Dygo Gdgg 

0.266 

0.22 

li^Frgg Dygg 

0.362 

0.31 

Wi06 ^^i°Hfi08 

46.2 

34.9 


by: 


two cases and the resonant state is an ex¬ 

cited state, (ii) In the other cases, the decay ^^^Gd is 
to a transitional nucleus, and the decays ^®"^Er and 
are to strongly deformed nuclei. For these reasons, one 
needs a model that can calculate reliably energies and 
wave functions of ground and excited states in spherical, 
transitional, and deformed nuclei. To this end we make 
use of the microscopic interacting boson model (IBM-2) 
(40| . The method is described in Refs. P, 01 . We write 


r<ECEC 

'^Ou 


1 


{GcosOY 


4R2 In 2 


27r^ 
2 


{h(?){meC 


, 2\7 


1 ’ ’ 


( 9 ) 


The obtained prefactors are listed in Table m 

Recently Krivoruchenko et al. [2l| have presented a 
theory of PF slightly different from ours. The PF ex¬ 
tracted from their paper is given in the last column of 
Table mi for comparison. Note the large value of 
in both calculations for decay, due to the large value 
of Za = 74/137 ^ 0.5. For these heavy nuclei the PF 
is very sensitive to the treatment of the electron wave 
function near the origin r = 0. Differences between the 
two PF calculations may in part arise from the way in 
which the nuclear size and electron screening correction 
is taken into account. 


III. NUCLEAR MATRIX ELEMENTS 

The calculation of NMF for the OvEGEG is more dif¬ 
ficult than for 0 jz/ 3“/3“ because of two reasons: (i) In 


with Mp‘^\ and defined as in Fq. (20) of 

Ref. @1 

= (^X;0+ \h^^[p)/g\fY-Jp) , 

= (^X;0+ \h^{p)/gl\^Y-JF) , (11) 

= (^X;0+ \h^{p)/g\\''Y-JF) , 

where h{p) are the form factors tabulated in Table II of 
01 . The wave functions of the initial and final states 
are taken either from the literature or from a fit to the 
observed energies and other properties {B{E2) values, 
quadrupole moments, B{M1) values, magnetic moments, 
etc.). The values of the parameters used in the calcula¬ 
tion are given in Appendix A. The quality of the wave 
functions as well as the quality of the description of ener¬ 
gies and electromagnetic transition rates is excellent both 
in transitional A=152, Fig. 01 and strongly deformed 
A=164, Fig. O and A=180, Fig. [51 nuclei. 

For decays to excited states, particular care must be 
taken in identifying the state in the calculated spectrum. 
For decay to ^®®Gd this state is the 0^ state calculated 
at 1988 keV, while for decay to ^^^Te this state is also 
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FIG. 4. Comparison between calculated and experimental low-lying spectra for ^®^Gd— 



FIG. 5. Comparison between calculated and experimental low-lying spectra for 




the 0^ state calculated at 2790 keV. The quality of the 
excited spectrum of states in ^^^Te, Fig. [3 and ^^®Gd, 
Fig. [ 8 ] is also excellent. 

Using IBM-2 wave functions and the theory previously 
described da we can calculate the NME for neutrino¬ 
less ECEC decay shown in Table IIIII The values of the 
nuclear matrix elements to the excited states are consid¬ 
erably smaller than those to the ground state, due to the 
very different nature of these states. To illustrate this 
point we show in Table m the results of the calculation 
for A = 156 leading to the first five O’*" states in ^®®Gd. 
We see that the matrix elements of 0;^, 0^, and 0^ are 
one order of magnitude smaller than Of and Of. This has 
a major consequence on QvECEC essentially excluding 
as possible candidates all decays to excited states. A 


similar conclusion was drawn in Ref. for the decay 

96Ru^.96Mo*. 

NME to the g.s of i52Sni, ^nd i 8 °Hf for light 

neutrino exchange have also been calculated in a vari¬ 
ety of methods. In Table El we compare our results to 
these other calculations. While the QRPA matrix ele¬ 
ments are of the same order of magnitude as IBM-2, the 
EDF results are much smaller, in particular they show a 
large reduction for the deformed nuclei ^®^Dy, and ^®°Hf. 
The origin of this discrepancy between IBM-2/QRPA 
and EDF is not clear. As shown in Fig. [B] the quality 
of the spectrum (and the electromagnetic transitions not 
shown) in IBM-2 is excellent including the location of the 
j3 and 7 bands, which is crucial for providing good wave 
functions of collective states. Also the collective struc- 
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FIG. 6. Comparison between calculated and experimental low-lying spectra for 
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FIG. 7. Comparison between calculated and experimental low-lying spectra for 




TABLE III. IBM-2 nuclear matrix elements (dimensionless) for neutrinoless ECEC decay with Argonne SRC and 

gv/QA = 1/1.269. 




light 




heavy 











i24xei24Te*[o+] 

i52Gdo-i52Sni[0)^] 

156j3y^l56Q^*[0j] 

i64Er^i64Dy[0+] 

18O.yy^l8Ojjf[0+] 

0.277 

2.132 

0.265 

3.456 

4.117 

-0.051 

-0.352 

-0.048 

-0.444 

-0.566 

-0.012 

0.095 

0.017 

0.221 

0.204 

0.297 

2.445 

0.311 

3.952 

4.672 

6.447 

68.910 

10.530 

107.900 

118.700 

-2.921 

-31.400 

-4.739 

-46.820 

-53.320 

-1.521 

12.810 

2.616 

32.870 

28.200 

6.740 

101.200 

16.090 

169.800 

170.900 


ture of the parent nucleus and daughter nucleus 


isoHf is practically identical, with a large overlap, and 
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FIG. 8. Comparison between calculated and experimental low-lying spectra for 


TABLE IV. IBM-2 nuclear matrix elements (dimensionless) for A = 156 neutrinoless ECEC decay to the first five O"*" 

states with Argonne SRC and gvjgA = 1/1.269 . 





light 




heavy 










Of 

2.796 

-0.398 

0.132 

3.175 

82.560 

-36.990 

17.480 

123.000 

02+ 

1.532 

-0.227 

0.076 

1.749 

47.630 

-21.360 

10.430 

71.320 

03+ 

0.403 

-0.065 

0.022 

0.466 

13.900 

-6.242 

3.225 

21.000 

of 

0.265 

-0.048 

0.017 

0.311 

10.530 

-4.739 

2.616 

16.090 

05+ 

0.302 

-0.046 

0.016 

0.346 

9.809 

-4.402 

2.204 

14.750 


TABLE V. Comparison between IBM-2 matrix elements with 
Argonne SRC for OvECEC decay and QRPA and EDF. 


Decay 



(light) 




spherical deformed 



IBM-2 

QRPA^ 

QRPA^ EDF^ 

^t'Gdss 

^ei^Smgo 

2.44 

7.59 

3.23-2.67 1.07-0.89 

^i'Erge 


3.95 

6.12 

2.64-1.79 0.64-0.50 

^fWi06 

^^i°Hfi08 

4.67 

5.79 

2.05-1.79 0.58-0.38 


Ref. m 
Ref. 0 


both nuclei have many collective pairs, =4 or 5 and 
Ni, =10 or 9, respectively, that contribute to the decay. 
Although there is a reduction in the NME due to the de¬ 
formation, we still would expect the matrix elements to 
be comparable to those in spherical nuclei. 


IV. HALF-LIVES 


The results for PF and NME of the previous sections 
can be combined to give the half-lives in Table IVII and 
Figs. [5] and (TUI The best case scenarios appear to be 
i80W^i®°Hf[0)^], for which, how¬ 
ever, the predicted half-life even with unquenched value 
of gA = 1.269 is the order of lO^^yr. 

As in the previous papers, an important question is the 
quenching of the gA which occurs to the fourth power in 
Eq. (2). To estimate this effect we use the parametriza- 
tion of Eq. (40) of Ref. [3|, (maximal quenching) 


^ i,269A-°-1®. (12) 


With this parametrization we obtain the values of Table 
EH (right). 
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FIG. 9. (Color online) Expected half-lives for light neutrino 
exchange with {rrijj) = 1 eV, qa = 1.269. The figure is in 
semilogarithmic scale. 



FIG. 10. (Color online) Expected half-lives for heavy neutrino 
exchange with irip = 2.75 x 10“^, qa = 1.269. The 

figure is in semilogarithmic scale. 


I80’v^^i80jjf[o+] and i52Gd^i52gj^[o+], where for light 
neutrino exchange ti /2 = 3.44 x lO^^yr and 8.03 xlO^^yr, 
respectively, for (m^) = 1 eV and gA = 1.269. For 
comparison decays have T-yj^ of the order of 

lO^^yr. The half-lives for OvECEC are, however of 
the same order of magnitude of the best case 

for this decay being that of which has 

Ti /2 = 3.32 X lO^^yr. Our conclusion is that, even in 
the optimistic scenario of gA = 1.269, OvECEC is un¬ 
reachable in the present generation of experiments. 

TABLE VI. Left: Calculated half-lives in IBM-2 for neutrino¬ 
less double electron capture with light neutrino exchange for 
{nil,) = 1 eV and heavy neutrino exchange with m,p (m~^^ = 
2.75 X 10“^, and using unquenched value qa = 1.269. Right: 
Same as right but using maximally quenched value g^A^ff^ ~ 
1.269A-“-i®. 




n/2(l0^^)yr 



unquenched 

maximally quenched 

Nucleus 

light 

heavy 

light 

heavy 

i24Xe^4Te*[0|[ 

1520 

150 

49000 

4820 


8.03 

0.237 

298 

8.83 

156Dy^l56Gd* [0|] 

2890 

54.7 

110000 

2080 

i64Er^i64Dy[o+[ 

1880 

51.6 

74000 

2030 

18O’VV^18OHf[0+] 

3.44 

0.131 

144 

5.47 
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VII. APPENDIX A: PARAMETERS OF THE 
IBM-2 HAMILTONIAN 


V. CONCLUSIONS 

In this paper, we have presented both PFs and NMEs 
for neutrinoless double electron capture and from these 
calculated the expected half-lives for both light and heavy 
neutrino exchange. The values obtained are longer than 
those for 0vP~f3~ decay. The best cases are those of 


A detailed description of the IBM-2 Hamiltonian is 
given in (H and Hg. For i24Xe, i52Gd, 156 q^^ 

and ^®^Gd the Hamiltonian parameters are taken from 
the literature. The new calculations are done using the 
program NPBOS [i^l adapted by J. Kotila. The values 
of the Hamiltonian parameters, as well as the references 
from which they are taken, are given in Table IVHI 
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